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^ . Abstract 

Motivated by a recent paper that compares the results of the analysis of the CKM 
angle a in the frequentist and in the Bayesian approaches, we have reconsidered the 
information on the hadronic amplitudes, which helps constraining the value of a 
in the Standard Model. We find that the Bayesian method gives consistent results 
irrespective of the parametrisation of the hadronic amplitudes and that the results of 
the frequentist and Bayesian approaches are equivalent when comparing meaningful 
probability ranges or confidence levels. We also find that from B — > tttt decays alone 
the 95% probability region for a is the interval [80°, 170°], well consistent with recent 
analyses of the unitarity triangle where, by using all the available experimental and 
theoretical information, one gets a = (93 ± 4)°. Last but not least, by using simple 
arguments on the hadronic matrix elements, we show that the unphysical region 
a ~ 0, present in several experimental analyses, can be eliminated. 



1 Introduction 



Motivated by the criticisms recently appeared in ref. [1], we present a new analysis of the 
CKM angle a from B — > nn decays based on the Bayesian statistical approach. The main 
results are the following: 

• we show that the differences between the frequentist and the Bayesian approaches 
are NOT due to the difference in the two methods but to the difference in the 
physical assumptions on the weak amplitudes, contrary to the claims of ref. [1]; 

• although we expect an eightfold ambiguity for a using as "a priori" knowledge 
only isospin symmetry and the experimental measurements of the relevant branch- 
ing fractions and CP asymmetries, this degeneracy can be reduced in the presence 
of further information on the hadronic amplitudes. We will discuss and use this 
information which can already be extracted from the data; 

• among the information that we do have on the amplitudes, the existence of a scale 
of strong interactions and the use of the experimental values of related decays allow 
to limit the size of the hadronic matrix elements and to eliminate the unphysical 
region of values corresponding to a close to zero; 

• within our approach, we obtain consistent results at a meaningful probability value 
irrespective of the parametrisation used for the hadronic amplitudes. 

Our analysis is performed within the Standard Model. All the technical details which 
are not relevant to the present discussion, but are used in the analysis, can be found in 
previous publications of the XJTfit Collaboration [2-5]. 

The paper is organised as follows. We first recall the basic formalism and definitions 
of the different quantities entering the analysis; we then discuss the case in which no 
further assumptions are made for these decays besides the experimental knowledge of the 
relevant branching fractions and CP asymmetries, and SU(2) isospin symmetry. In this 
case we make a comparison between frequentist and Bayesian methods along the lines of 
ref. [1]; we finally discuss the present experimental information on the weak amplitudes, 
and the corresponding constraints, and show how this information helps in reducing the 
ambiguity of the solution. 

We think that neither our study nor the study of ref. [1] can be decisive in settling 
the long standing struggle on the validity of the frequentist and Bayesian methods. For 
this reason, contrary to what was done in ref. [1], we do not find very useful to present 
in a physics paper a long list of learned citations, with philosophical statements and 
sentences taken from illustrious references in favour of the Bayesian statistics. Rather, 
we thank the authors of ref. [1] who stimulated us to reconsider the analysis of the CKM 
angle a, and allowed us to improve our analysis by using the available experimental 
information. The important point is not which of the two approaches is used rather the 
ability to get predictions which can eventually be verified "a posteriori". The success 
of a phenomeno logical analysis relies on the capacity of anticipating the correct result. 
Among the successful (and rather accurate) predictions of our collaboration let us recall 
sin 2(3 [2,6] and ArriB 3 [2,3,7], made long before their experimental measurements. 
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SET 1 ref. [1] 


5.5 ±0.6 


5.1 ±0.4 


1.45 ±0.29 


-0.37 ±0.10 


-0.50 ±0.12 


-0.28 ±0.40 


SET 2 ref. [8] 


5.7 ±0.4 


5.2 ±0.2 


1.31 ±0.21 


-0.39 ±0.07 


-0.59 ±0.09 


-0.37 ±0.32 



Table 1: Experimental numbers used in the analysis of the CKM angle a. The branching 
fractions are given in units of 10~ 6 . 

2 Generalities 

The physical quantities which are measured in the experiments and used in the extraction 
of a from B —> ixix decays are the CP-averaged branching fractions £>+°, B+~ and B®° n and 
the coefficients , and S+~of the time-dependent CP asymmetries. At present <S°° 
is not measured. The branching fractions and CP asymmetry coefficients are labelled by 
the charge of the final state pions. We give in Tabled] the two sets of values that we have 
used in this analysis: SET 1 is the set of ref. [1], which will be used for a comparison of 
the results; SET 2 is taken from HFAG [8] and contains the values updated at the summer 
2006 conferences. The CP asymmetry coefficients and branching ratios are related to the 
amplitudes by the relations 

C% = (\A ij \ 2 - \A ij \ 2 )/(\A ij \ 2 + \A^\ 2 ) , SH = -2 lm[A ij A ij *]/(\A ij \ 2 + |^'| 2 ) , 
^V°° = (l^'T + |A + ->°T)/2 , B£ = t b+ /t b o (\A +0 \ 2 + \A-°\ 2 )/2 . (1) 

In the equations above we have included in the definition of the amplitudes trivial and 
standard factors such as the two body phase space and the squared Fermi constant. Here 
and in the following the branching fractions are given in units of 10~ 6 . In our convention, 
denoted as "natural units" in the following, the factorised amplitudes are of 0(1), as 
explained below. 

In the Standard Model, a minimal but non zero set of "a priori" knowledge on strong 
interactions is common to all the phenomenological analyses of a. The universal "a priori" 
assumptions of these studies are that strong interactions are flavour independent and 
conserve parity and CP. This information is born out from the experimental measurements 
of strong interaction processes. In addition most of the analyses are performed in the 
approximation in which isospin symmetry breaking effects - including electromagnetic 
corrections - are neglected. □ In the remaining of this paper we will denote this set of 
"a priori" assumptions, combined with the isospin symmetry approximation, "minimal 
assumptions" . 

Using the "minimal assumptions" , in the standard parametrisation the amplitudes are 
written as [9] 

A + - = A(B° -> 7T+7T-) = e ~ ia T + ~ + P, 

A 00 = A(B° - vrV) = -L ( e --T 00 - P) , (2) 



A +u = A(B + -> tt+tt ) = ^~ ia (T 00 + T~ 



1 For lifetimes we will use the experimental value for the charged and neutral B mesons 
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The six free parameters are the absolute values of P u and P, \T^\ and |P|, their rel- 
ative strong phases, <pp =Arg[T + ~ P*\ and 0o =Arg[T +_ T 00 *], the overall phase being 
irrelevant, and a. 

Starting from the general expressions in Equation (J2J), under the "minimal assump- 
tions" mentioned above, one finds for a [9] either zero or eight solutions from Equa- 
tions ([T]), corresponding to [1,10] 
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We also give the formulae from which it is possible to extract, from the branching fractions 
and the CP asymmetry coefficients, the hadronic parameters T % \ P and the relative phases 
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arctan(xj, is the value of the angle fa obtained when both yi 



T€m.d>i and 



Xi = cos (f>i are known (with T a positive factor). In the previous equation for |T + ~ 
and |T 00 | one must choose the unique combination of signs that reproduces the correct 
result for B +0 

It has been argued in ref. [1] that using the "minimal assumptions" defined above, the 
frequentist and Bayesian approaches give different results. We show that this is not the 
case with two clear examples, one based on the present experimental information, and a 
second one in the hypothesis that the experimental errors are reduced by a factor ten. For 
this comparison we use the same data as in ref. [1], SET 1 of Table [TJ In Figured] we show 
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the frequentist and Bayesian results with the present experimental errors. To obtain 
these plots, in the Bayesian approach branching ratios and CP asymmetry coefficients 
are extracted with flat a priori probability density function (p.d.f.) and weighted by the 
experimental likelihood. The ES parametrisation defined in ref. [1] is the only one where 
no priors for the hadronic parameters are specified. Therefore this is the only case in 
which no additional physical information with respect to the frequentist fit is used. 

Very often there is some confusion in the interpretation of the results because the 
shape of the frequentist and Bayesian figures look very different. It is important to stress 
that this happens because they correspond to two different quantities. Therefore it may 
be useful to recall how to read the relevant physical information in the two cases. Indeed 
the shape is not the important issue to this purpose. In the Bayesian case what it is 
usually shown is the p.d.f. for a given physical quantity, whereas in the frequentist case 
what it is shown is the confidence level (C.L.), or rather 1-C.L., corresponding to a certain 
value of the physical quantity. Although there is not a rigorous correspondence, in order 
to compare the two approaches one might confront the 68% (95%) integrated probability 
region of the Bayesian case with the set of values corresponding to the 68% (95%) C.L|j 
In the case of Figured] we have explicitly indicated with a dark band the 95% integrated 
probability region and the region of values corresponding to the 95% C.L. in the Bayesian 
and frequentist case respectively. The same considerations apply to the figures shown in 
the remaining of this paper. 

From the selected regions in Figure[T]our conclusion is indeed just the opposite of what 
is claimed in [1]: at a confidence level which is meaningful for obtaining some physical 
information, for example at 95% C.L., and even at the 68% C.L., when using the same 
physical assumptions (in the present case the "minimal assumptions"), the frequentist 
result is equivalent to the Bayesian expectation, namely the range of values of a between 
25° and 65° is excluded whereas the complementary interval is fully allowed. Even in the 
frequentist case, the eight solutions emerge only for values of the C.L. smaller than 5%, 
which have no physical meaning. This is further illustrated by studying the case in which, 
at fixed central values, the errors are reduced by a factor of 10: as shown in Figure [2] the 
eight solutions are separated both in the Bayesian and frequentist case. Note that even 
in the frequentist case they are still grouped at 95% C.L.. Thus we conclude that there 
is no substantial difference in the physical information obtained in two approaches, as 
also discussed in many other examples and stated in the Yellow Report of the first CKM 
workshop [11]. 

To strengthen even more the previous arguments we have repeated the frequentist fit 
with the updated values corresponding to SET 2 of Table dj taken from HFAG [8]. The 
results are shown in Figure [3J In this case, the new values of the branching ratios make 
the eight solutions overlap also in the frequentist approach (independently of the chosen 
confidence level). This shows that the separation of the eight solutions in the frequentist 
approach was just a fortuitous accident which may disappear for small changes in the 
central values, when the errors are large. 

We have also repeated the exercise of removing some crucial experimental information, 
namely we ignored the experimental measurement of [1]. In the Bayesian approach 
this corresponds to let £>°° free to vary between zero and 10 6 . In this case we get for the 

2 Similar figures were shown in Figure 2 (first and fifth) of ref. [1]. 

3 In our analysis, probability regions are delimited by the intercept of the distribution with a horizontal 
line. 
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Figure 1: Comparison of Bayesian (left) and frequentist (right) result in the "minimal 
assumptions" case. The dark band in the left figure corresponds to the 95% probability 
region, to be compared with the 95% C.L. interval obtained in the right figure. The same 
interval is essentially selected in the two cases. The same conclusion holds when comparing 
the 68% probability region to the 68% C.L. interval. 

p.d.f. of a the result shown in Figure H] which allows, for the 95% total probability region, 
almost all the values between 0° and 180° (corresponding to the 95% of the total area 
under the curve), fully consistent with the frequentist approach. We could not reproduce 
Figure 4 of ref. [1] (case ES) which we think contains some error producing the deep at 
a ~ 45°. 

For the remaining of this paper, in the analysis we use the more recent values of SET 

2 in Table [H taken from HFAG [8]. 

3 Information on Hadronic Matrix Elements 

In this Section we present several arguments which show that the size of the hadronic 
matrix elements of the operators appearing in the effective weak Hamiltonian is indeed 
of the order of magnitude which is expected in QCD. This implies that T % i and P are 
numbers of 0(1) in the natural units used in this paper and that the range of values 
we used in our previous analyses was too pessimistic and can be restricted. With this 
improvement, as shown in the next Section, the dependence on the prior noticed in ref. [1] 
is substantially eliminated and the constraint on a improved. 

As mentioned before, it is illusory to allege to be able to perform an analysis of 
a without a minimal "a priori" knowledge about strong interactions and the hadronic 
parameters. Such knowledge is encoded in the expressions given in Equation (j2J) which 
are valid only if strong interactions are flavour blind and CP conserving, besides using 
the approximation that isospin breaking effects are negligible, which implies Aqqd ^ 
(m d - m u ).0 

4 Electroweak penguin contributions and electromagnetic corrections to the decay amplitudes are also 
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Figure 2: Same as Figure U\ but with the experimental errors reduced by a factor 10. 
Obviously, also in the Bayesian case, an eightfold ambiguity appears. 

This however does not exhaust the information that we have on the size of the matrix 
elements. A general consideration is that we believe that the theory of strong interactions 
is QCD, which is a renormalisable theory with a dimensionless coupling constant, and has 
a natural scale of 0(1 GeV). Thus, on dimensional grounds, in the absence of other scales 
we expect {MiM2\0\M) ~ A.q CD for a dimension-six four fermion operator oEI We would 
be very surprised to find that this matrix element has a size of the order of (ITeV) 3 or 
M pianck- m the presence of other scales, such as the mass of a heavy quark, the situation 
is slightly more complicated, but the argument is conceptually similar. Indeed we expect 
that [13]- [15] 

/A \ 3 / 2 

(mc\d\B) ~ UM B f+(0) ~ f n M B (-Jg£j ~ M B /2 \f CD . (5) 

Note that this simple dimensional argument, which we have derived using f actor ised ex- 
pressions in the intermediate steps, has a more general validity than factorisation al- 
though, obviously, factorisation respects the same scaling laws. For chirally enhanced 
contributions, the enhancement factor, by which the expression in Equation (jSJ) must be 
multiplied, is of the order 2M^/(m u + m ( i)/MB ~ 0.8 and thus it does not change the nat- 
ural size of the amplitude. It is straightforward to show, see Equation ([6]) below, that this 
implies T u ~ 1. Again we would be very surprised to find that, in order to reproduce the 
experimental branching ratios, the size of the matrix elements must be much larger than 
its natural size. In the language of factorisation this would correspond to dimensionless 
i?-parameters much larger than one. 

We now provide further support to these general considerations: 

1. The first exercise, useful to estimate the range of values which can be expected, is to 

ignored. See ref. [12] for a recent discussion of isospin breaking effects in the extraction of a. 
5 For example (bu)v-A (ud)v-A hi the notation of refs. [16]. 
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ref. [17] 


ref. [18] 


ref. [19] 


Exp. 


3.6-5.3 


4.3(1 ±0.3) 


3.711:? 


5.5 ±0.6 



Table 2: Comparison of the experimental value of B^® with the theoretical predictions 
[17-19]. All values are given in units o/10 -6 . 

compute T + ~ using only current- current operators and strict factorisation, namely 
|T+-| 2 = ^^^. \C x {M B ){^-\O x \E^ + C 2 (M B )(n+n-\0 2 \B°)f x 10 6 



?2 ^_|T/ t/»|2 ^ 2 



uii I 



32?rM B 



x|MiA/ + (0)| 2 xl0 6 . (6) 



Using f n = 132 MeV, /+(0) = 0.3, Ci(M B ) = -0.2 and C 2 (M fl ) = 1.1, we find 
|T + ~| = 3.2 in natural units. 

2. Long before it was experimentally measured, several predictions based on factorisa- 
tion or other approaches existed for , £>°° and £>+°. The latter is the simplest 
to predict since it only depends on emission diagrams, without penguins, annihila- 
tions or further complications. In Table [2] we compare some predictions, obtained 
in different theoretical frameworks, with the present experimental determination. 
We find that the order of magnitude of the predicted branching fractions, based on 
values of the hadronic amplitudes of 0(1) in natural units, are approximately in 
agreement with the experimental value. 
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Figure 4: Probability distribution for a obtained removing . No significant information 
can be obtained in this case. 



3. The third argument in favour of the natural order of magnitude for the matrix 
elements is based on the scaling of the rates between B and D decays. Since, in the 
heavy quark limit, the partial decay rates scale as the squared amplitude divided 
the heavy mass, namely as 1/M x \M 1 / 2 KqI ;d \ 2 ~ ^qcd-> ^ we take the ratio 

|r+-(j^7r+7r-)| 2 \V ub V: d \ 2 

iT+-p°-rt-)i 2 ~ ivw 1 1 

this is independent of the meson mass and can be used to extract the absolute value 
of the unknown hadronic parameter T H as follows 

\T+-\ 2 = BR(D° -> 7T+7T-) x 10 6 ^-i?. (8) 

Using the central value of the experimental measurement BR(D° — > 7r + 7r~) = 1.5 x 
1(T 3 , t d = 0.41 x 1(T 12 sec, r B = 1.6 x 1(T 12 sec, \V ub \ = 3.7 x 1(T 3 , |1/ cd | = 0.22, 
we find |T H | = 1.3. 

4. We can get some knowledge about the parameter P from the study of B s — > K + K~ 
decay. Up to doubly Cabibbo suppressed terms, this decay proceeds only through 
the penguin contribution P s , which corresponds to P up to SU(3) breaking ef- 
fects [20]. Using the central value of the experimental measurement BR(B S — > 
K+K-) = (24.4 ± 1.4 ± 4.6) x 10~ 6 [21], and the relation 

\P\ 2 = BR(B S - K + K~) x 10 6 > ( 9 ) 

1~B° \VtsV tb \ 2 
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we find \P S \ = 1.1. Even accepting the very pessimistic point of view that SU{3) 
breaking effects are 100%, one still obtains a number of 0(1) which automatically 
constrains, when combined with B — > tttt decays, also T H to be of (9(1). For B s — > 
K + K~ this argument could be spoilt if the Cabibbo suppressed emission amplitude, 
T, were much larger than the penguin amplitude, namely if \T\ ^> \P S \ ^ 1 in 
natural units. This possibility is however excluded because, in order to reproduce 
the experimental values of £>+~ and B®° n , we may have \T\ 1 and |P| ^> 1 but it 
is necessary that |T| ~ \P\. 

In summary all the experimental and theoretical evidence is in favour of matrix elements 
of 0(1) in "natural units". The arguments from 1. to 3. show that in our previous 
analyses [4,5] the range that we used by varying |T lJ '| and \P\ between and 10 was 
rather conservative and easily accommodates 0(1) corrections to the theoretical estimates 
presented above. Argument 4. puts a stronger constraint on the penguin contribution 
since, even assuming an SU(3) breaking effect as large as 100%, it is difficult to imagine 
that it could exceed the value \P\ = 2.5 in natural units. As discussed in 4. this upper 
value on \P\ limits the acceptable values of |T +-,00 |. 

We conclude this Section with a few comments on appendix B of ref. [1]. In this 
appendix, the authors criticise the Bayesian approach for being unable to reproduce the 
peak of the p.d.f. at a around zero. We believe instead that they are completely mislead 
by their prejudice of ignoring the existing information on the Standard Model and on the 
hadronic matrix elements and for this reasons they are unable to eliminate the unphysical 
solutions at a close to zero, which cannot be there. A very simple argument kills 
simultaneously the solutions at a ~ and their arguments to show that the Bayesian 
approach is unable to cope with the RI parametrisation, where T + ~ = r + ~/ sin a and the 
prior distribution is taken for r +_ rather than for T + ~ . For example, it is straightforward 
to show that a < 2° corresponds to T + ~ > 30 in natural units. This implies that, in order 
to fit the branching ratios, also P ~ 30, i.e. an SU(3) breaking effect of about 3000%! 
Bearing in mind that m s /md ~ 10 we would then expect SU{2) effects of "only" 300% 
which should, in our subjective and humble opinion, be taken into account invalidating 
the assumption of neglecting isospin breaking effects in the analysis. 

4 Adding Useful Information 

In Section [2] we have shown that, using the "minimal assumptions" on the hadronic 
matrix elements, no real difference in the physical information at a significant level of 
confidence exists between the Bayesian and the frequentist approach; in Section [3] we 
have discussed with several examples the information that we have on the weak hadronic 
matrix elements. All the arguments support the existence of a typical size for the hadronic 
matrix elements, which constrains the range of possible values. We do not understand 
why one should ignore this knowledge while using the simplification (and the information) 
that comes from isospin symmetry. 

In this Section, we make use of the constraints on the size of the matrix elements in 
different ways. In general, we show that the information on the matrix elements helps 

furthermore they find the rather original result that the C.L. is about 90% for a ~ 1CP 30 and zero 
for a = 0!! 
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eliminating some of the eight solutions that exist in the "minimal assumptions" case. The 
prejudice that any acceptable method must lead to eight solutions [1] is obviously wrong: 
if we could determine exactly the absolute values of the hadronic matrix elements and the 
strong phases from a theoretical calculation, from the experimental values of £>^4 and of 
the CP asymmetry coefficients we could extract unambiguously the value of a. A partial 
knowledge, or some constraint, will then remove at least partially the degeneracy of the 
solutions. 

In particular, we show that, using the constraints on the hadronic matrix elements, 
not only we eliminate the pathological solution at a ~ 0, which could survive only with 
stratospheric values of the matrix elements, but we can substantially reduce the differ- 
ence between the different parametrisations that was emphasised in ref. [1]. Indeed the 
difference disappears for the 95% probability region, and is marginal even for the 68% 
one. 

For the parametrisations we use the same convention as in ref. [1]: MA in which the 
absolute values of the hadronic amplitudes and their phases are extracted with a flat p.d.f.; 
RI in which the real and imaginary part of the hadronic matrix elements are extracted 
with a flat p.d.f.; ES in which, in the Bayesian case, we generate the measured quantities, 
namely £>+°, £>+~ and £>°° and the CP asymmetry coefficients C£~, and S£~, with 
flat a priori probability density function (p.d.f.) and weighted with the experimental like- 
lihood. In this case, when solving for the absolute value of the matrix elements according 
to Equations (jlj) of Section [2j we only accept the solutions when they fall in the "a priori" 
acceptable region defined below. 

For the following discussion, we distinguish the "previous" set of a priori, which cor- 
responds to the ranges used in our previous analyses (|T y '| < 10, \P\ < 10, and arbitrary 
phases) [4] and the "current" one, corresponding to \T l ^\ < 10, \P\ < 2.5 and arbitrary 
phases. In the last case we limit \P\ to be less than 2.5 since, to our knowledge, no physical 
quantity suffers from a SU(3) breaking effect larger than 100% and from B® K + K~ we 
estimated \P\ ~ 1.1. In this Section we always use the values of the measured quantities 
corresponding to SET 2 in Table [TJ 




o 

I— 

a. 





Figure 5: Probability distribution for a obtained with three different prior assumptions: ES 
(left), MA (centre) and RI (right). In all the cases the hadronic parameter are constrained 
by |T" U | < 10, \P\ < 2.5. In all the cases the region with "low" a is excluded. This is 
not due to the statistical approach, but to the physical assumptions on the hadronic matrix 
elements (the smooth aspect of these p.d.f.s is due to the fact we use only 10 5 (sic!) Monte 
Carlo events supplemented with an efficient numerical algorithm). 
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Figure 6: We show the range of allowed values of a corresponding to the 95% probability 
regions for the three parametrizations (MA=red, RI=green, ES=blue) as a function of the 
maximum allowed value of \P\, Pmax, scanning the range 1 < Pmax < 5.5. This figure 
shows that the residual dependence on \P\ is very mild. 

In Figure [5] we show the probability distribution functions for the three parametri- 
sations (MA, RI and ES), in the "current" case. A comparison of the 95% probability 
region, highlighted in dark, show that the different parametrisations give the same physi- 
cal information, although the shape of the p.d.f.s looks different. The difference observed 
in ref. [1] between MA and RI in the "previous" case was due to the different "a priori" 
distribution for the absolute values of the amplitudes (RI has a p.d.f. linearly grow- 
ing with the absolute value) and the lack of further information on the hadronic matrix 
elements. In this respect we have to admit that we were not Bayesian enough in our previ- 
ous analyses, because the "previous" upper limit consisted in taking the "intuitive" upper 
bound one order of magnitude larger than the typical size, without using the experimental 
information which was already available. 

Indeed the upper bound on P, which was obtained using only SU(3) arguments, 
implies a bound also on |T lJ '|. If P is limited, in order to reproduce the observed B — > tttt 
branching ratios, |T U | cannot be too large. We leave to our friends of ref. [1] the exercise 
to show that with the same physical inputs one gets in the frequentist case the same 
results at the 95% C.L.. 

After the above discussion, where we have shown that the Bayesian approach gives 
physically meaningful results, we are now ready to discuss the implications of the bound 
that we have imposed on the size of the hadronic amplitudes, in particular the penguin 
amplitude \P\. 

Even in the case of observed CP violation in B — > tttt decays, the system of equa- 
tions used to extract a and the hadronic amplitudes, including the strong phases, from 
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branching fractions and CP asymmetry coefficients, Equations ([3]) and (J4|), admits the 
unphysical solution a = 0. This solution can only be obtained however in the peculiar 
limit \T\ — > oo, \P\ — > oo and P/T — > — 1. This has several unappealing features: i) in 
the Standard Model, CP violation must disappear for a — > since the unitarity triangle 
collapses to a line and the Jarlskog determinant vanishes; ii) the hadronic amplitudes 
cannot go to infinity without violating the basic properties of any renormalisable field 
theory (and we have discussed what is their expected size); iii) the unphysical solutions 
for a ~ are eliminated, since we have shown that they imply arbitrary large ST/ (3) 
breaking effects, which are not consistent with the assumption of isospin symmetry in the 
expressions of the amplitudes. 

The bounds on the hadronic amplitudes discussed in Section [3] allow to get rid of the 
unphysical solution in a straightforward way. Rather than living with a non physical 
solution, we prefer by far to remain with the smallish residual dependence on the upper 
value on \P\ on the allowed interval for a, shown in Figure [61 It would be interesting 
to compare the uncertainty originating from the upper value on \P\ to the size of the 
uncertainty due to the corrections coming from the neglected isospin breaking effects. 

For completeness, we finally give in Figure[7]the p.d.f.s for the relevant amplitudes and 
strong phases, as selected by our analysis, in the MA parametrisation. These distributions 
can be used for comparison and test of the different theoretical approaches which compute, 
within some approximation, the hadronic matrix elements from QCD [18,19,22-24]. We 
present two cases: the case in which only the quantities measured in B — > tttt decays are 
used, and the case where, in the Standard Model, we make a combined CKM analysis using 
all the available experimental and theoretical information [2-5]. The results confirm the 
a priori knowledge on the weak hadronic matrix elements based on dimensional analysis 
and the existence of a typical scale in QCD, namely the allowed values of |T + ~ ,00 | and |P| 
are of 0(1). To be more precise, in the full Standard Model analysis, from the 68% region, 
we find |P| = 0.80 ± 0.24, |T + ~| = 2.1 ± 0.1 and |T 00 | = 1.4 ± 0.2, in full agreement with 
the expectations based on simple physical arguments of QCD and discussed in Section O 

Conclusions 

Stimulated by a recent paper on the extraction of the CKM angle a from B — > ttti de- 
cays [1], we have upgraded our Bayesian analysis with the following results: i) we have 
shown that the present information on the hadronic matrix elements, obtained from gen- 
eral theoretical arguments and experimental measurements, already allows a substantial 
reduction of the eightfold ambiguity in the determination of a, in particular by eliminat- 
ing the solutions at a ~ 0, that correspond to unphysical values of the amplitudes; ii) 
the information on the hadronic matrix elements substantially eliminates, in the Bayesian 
approach, the dependence of the results on the "a priori" probability distributions of the 
hadronic matrix elements, which was noticed in ref. [1]. Contrary to the claims of ref. [1], 
we have also shown that the differences between the frequentist and the Bayesian ap- 
proaches are not due to the difference in the two methods but to the difference in the 
physical assumptions on the weak amplitudes. We believe that a continuation of this ster- 

7 In addition to the arguments given in the text, the limit mentioned above would require an infinite 
amount of fine tuning and correlation among parameters related to different physics, namely a which has 
an electroweak interaction origin and the hadronic amplitudes which are sensitive to strong interactions. 
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Figure 7: We show the p.d.f.s for \P\, |T + ~ 
in the MA parametrisation, using only B - 
analysis (last two rows). 



and \T 00 \, and the relative phases, obtained 
mi decays (first two rows) or the full UT 



ile polemics in favour of the frequentist or Bayesian approach is only a waste of time and 
energies and it is better to concentrate the efforts in trying to combine in the most efficient 
way the rich information which is coming from the measurements of several non-leptonic 
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channels to constrain the CKM parameters and to make accurate predictions. 
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